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THE GEOPOTENTIAL AT SYNCHRONOUS-ORBIT ALTITUDES 

C. A. Wagner 


I would like to discuss briefly the latest solution I have obtained for the 
Earth’s gravity potential at synchronous-orbit altitudes. 

As you may know, a distant satellite in a synchronous orbit is affected 
significantly by only a few low degree terms in the conventional spherical 
harmonic representation of the field. These terms include the oblateness of 
the Earth and, because of the commensurability of the synchronous orbit, 
all the longitude-dependent terms as well. 

It has been 7 years since the first such satellite was put into orbit. Since 
that time, about 50 others have followed, making synchronous orbits 
probably the most popular variety of orbit. However, only about 15 of 
these objects are close enough to exact commensurability to be really useful 
in studies of the small but persistent effects of the so-called resonant 
harmonics of gravity. 

I have obtained the present solution for these harmonics from the most 
representative sample of tracking data from eight of these satellites in 21 
distinct free-drift arcs. The length of time for these arcs ranges from 6 weeks 
to 3 x /£years (Table 1). 

The observations I use in my solution are the mean Kepler elements for 
the satellites reported by the various responsible tracking agencies, including 
GSFC. The element most sensitive to the resonant effect is the mean 
geographic longitude or equator crossing of the satellite. This is a simple 
combination of the Kepler orbit arguments and the position of the Earth 
and is the principal observable recovered in the solutions. The quality of this 
observable varies from arc to arc, but in no case does it approach the fine 
precision of the tracking required for geodetic recovery from nonresonant 
satellites. The geodetic solutions with these data demonstrate conclusively 
that we are at last getting absolute determinations of the low-degree field. 
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This can be seen from the answers to three questions concerning these 
sensitive resonant data. First, if there were no resonant gravity harmonics, 
what kind of recovery of these data would we get? Figure 1 shows the 
evolution of the longitude of Syncom 2 between 1965 and 1969. Without 
resonant harmonics, this evolution would be almost linear, and, since the 
data in this arc have an accuracy of about 0.05 deg, the mean weighted 
residual in a best-fit trajectory would be over 100. Second, how well can the 
data be recovered with the use of recent gravity models that include 
resonant harmonics but not the data itself? Table 2 gives results of orbits for 
these satellites calculated with recent fields through fourth degree. The first 
two rows show the overall recovery with two recent Smithsonian geo- 
potentials. From about 100 without low degree coefficients, the mean 
weighted residual (rms) comes down to below 3 with the latest SAO field, 
which includes ground gravity information. An improvement in the fit to 
less than 2 is found with a 1969 all-satellite SAO field, very little different 
in these critical harmonics from the 1970 field. 

Third, what field adjustment is necessary for complete recovery of the 
data? This question is answered in the last two rows of Table 2. In the 
unconstrained solution, the dominant harmonics of second and third degree 
were freely derived along with the orbits. The residuals were lowest, but the 
(3,1) harmonic was not realistic and many correlation coefficients were 
large. The last row shows the most satisfactory solution obtained from the 
24-hour data. Here, all five significant harmonics were derived with 
reasonable a priori constraints from the SAO COSPAR coefficients. This 
constrained field shows adequate data recovery with realistic coefficients 
that are only moderately correlated. 

There should be little doubt that we can now predict the long term 
drift of these satellites to considerable accuracy even over periods of many 
months. In fact, a dramatic test of this proposition came in September 1970 
when Syncom 3 was reacquired by GSFC within 0.4 deg of its predicted 
position after being silent for almost 2 years. This prediction was made with 
the constrained field. 

The economics of computing these fields from mean element data is 
also interesting. Only 279 observations were used, compared with over 
50,000 for the SAO solutions. The solutions themselves involved about five 
differential corrections to both orbits and field parameters or a total of 60 
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orbit years of calculation, and took less than 4 min of computing time on an 
IBM 360/95 machine. 

The overall results of this study of the geopotential for synchronous 
satellites are as follows: 

1. Absolute accuracy of second-degree coefficients has been established 
to better than 3 percent. Coefficients through fourth degree have been 
established to better than 15 percent. 

2. The positions of equilibrium points for geostationary satellites have 
been established to better than H deg. 

3. Accuracy in prediction of orbits for 24-hour satellites has been 
established to better than 1 deg for periods greater than 2 years. 



Figure 1— Equator crossings for SYNCOM 2, 1965-1968. 
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Table 1 -Twenty-four hour satellite arcs in 1970 resonant geopotential 

solutions. 


SATELLITE 

ARC 

NUMBER 
OF KEPLER 
ELEMENT 
SETS 
USED 

SPAN OF 
DATA 
(modifitd 
Jut IX) 
0*yl) 

GEOGRAPHIC 
LONGITUDE 
SPAN IN 
ARC 

Idtfml 

ORBIT 

INCLINATION 

ESTIMATE OF 
0UALITY OF 
DETERMINATION 
OF THE GEOGRAPHIC 
LONGITUDE, n. {». 

IN ARC. FROM 
INDIVIDUAL ARC 
ANALYSES 
(ikgrMtl 

RMS RESIDUALS IN 
LONGITUDE FROM 
COMBINED ARC SOLUTIONS 

WITH 24 

HOUR SATELLITES 
ALONE. 

SOLVING FREELY 
FOR 6 RESONANT 
COEFFICIENTS 

WITH 24 HOUR 
SATELLITES 
COMBINED WITH 
SAO 1969 COSPAR 
DATA, IN CONSTRAINED 
SOLUTION FOR 10 
RESONANT COEFFICIENTS 

WITH SAO 
1969 COSPAR 
FIELD 
ALONE 

SYNCOM 2 

1 

IE 

36263 31351 

302 305 

33.0 

6 025 

0013 

0.013 

0.913 

SYNCOM 2 

2 

11 

313*1 38443 

296301 

32.8 

0 025 

0034 

0 034 

0.035 

SYNCOM 2 

4 

9 

36610- 34570 

196 243 

32.5 

0.020 

0.015 

0.014 

0.014 

SYNCOM 2 

S 

24 

30000-10007 

72109 

32.3 

0.635 

in 

0.031 

0.036 

SYNCOM 2 

1 

IS 

34416 34914 

65-60 

31.1 

0.020 

0.010 

0.016 

0.011 

SYNCOM 2 

DoO 

42 

30010-40104 

65 46 

79-32.0 

0.040 

0.037 

0.031 

0.103 

SYNCOM 3 

0 

9 

30699 36750 

170-140 

0.0 

0015 

0014 

0.015 

0.013 

SYNCOM 3 

7 

10 

36775 34035 

174-101 

00 

0.055 

0054 

0.065 

0.052 

SYNCOM 3 

11 

11 

39075-39262 

165172 

0.5 

0020 

0.012 

0.014 

0.013 

SYNCOM 3 

13 

4 

39345 39491 

160161 

1.3 

0.615 

0.006 

0.006 

0.010 

SYNCOM 3 

14 

IS 

39663-40175 

146160 

23.0 

0 045 

0.061 

0.061 

0.110 

ATS 1 

1 

12 

40204-40248 

210211 

1,2 

0.065 

0 001 

0.001 

0.001 

ATS 3 

1 

12 

4019740241 

200 312 

06 

0.005 

0002 

DM3 

0.003 

ATS 3 

2 

13 

4076740337 

717 700 

03 

B.M5 

0.004 

0.003 

0.002 

ATS 3 

3 

10 

40524-40577 

313-315 

0.3 

0.005 

0.002 

0.003 

0.003 

INTELSAT 2-F3 

1 

12 

39607 39905 

349 352 

1.0 

0 025 

0026 

0.021 

0.031 

INTELSAT 2 F3 

2 

10 

4040640642 

346349 

1.0 

0 020 

0.011 

0.012 

1.015 

INTELSAT 2 F4 

1 

1 

40323 40004 

179-194 

10 

0040 

0.043 

0.046 

• M2 

EARLYBIRD 

1 

11 

30497-39000 

324 332 

0.5 

0 025 

0.025 

0.026 

0.026 

EARLYBIRD 

2 

6 

39096-39219 

321-331 

0.1 

0.010 

0 007 

0.009 

0.009 

ATS 5 

1 

14 

4047640550 

253 256 

26 

0005 

0.004 

0.006 

9.009 


Table 2-Resonant geopotential coefficients for 24 hour satellites and rms 

fits to 21 arcs. 


FIELD 

OVERALL 
WEIGHTED 
RMS IN 
21 ARC 
SOLUTION 

2.2 

C 

S 

3.1 

C 

S 

3,3 

C 

S 

C 

4,2 

S 

4,4 

C S 

STANDARD 
SAO EARTH H 
(1970) 

2.46 

1558 

0 881 

2.128 

0.281 

0096 

0.199 

0.074 

0 158 

-0.0017 0.0072 

SAO 

COSPAR 

(1969) 

1.38 

1 566 - 

0.896 

2.040 

I 

0.262 

0.096 

0.198 

0.073 

0.148 

-0.0028 0.0078 

WAGNER 

(UNCONSTRAINED) 
2,2-3, 1 3,3 
(1970) 

0 849 

1 568 - 

■0.907 

! 

1.687 

0.483 

0 103 

0.204 

0.074 

0 158 

-0.0017 0.0072 

WAGNER SAO 
COSPAR 
(CONSTRAINED) 
(1970) 

0 896 

1.570 - 

■0.908 

2029 

0.289 

0 098 

0.205 

0.075 

0 150 

-0.0028 0.0078 


UNNORMALIZED COEFFICIENTS IN UNITS OF 10' 6 





